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K.^ ' We report a detailed theoretical investigation on electrochemical capac- 
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t:;;j- ' itance of a nanoscale capacitor where there is a DC coupling between the 

cn ■ 

^^ ' two conductors. For this "leaky" quantum capacitor, we have derived general 



analytic expressions of the linear and second order nonlinear electrochemical 
C^ ' capacitance within a first principles quantum theory in the discrete poten- 

tial approximation. Linear and nonlinear capacitance coefficients are also 
derived in a self-consistent manner without the latter approximation and the 
self-consistent analysis is suitable for numerical calculations. At linear order, 



^ ' the full quantum formula improves the semiclassical analysis in the tunnel- 



ing regime. At nonlinear order which has not been studied before for leaky 

capacitors, the nonlinear capacitance and nonlinear nonequilibrium charge 

show interesting behavior. Our theory allows the investigation of crossover of 

capacitance from a full quantum to classical regimes as the distance between 

the two conductors is changed. 
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I. INTRODUCTION 

The most significant development in electronic devices has been the progressive miniatur- 
ization of them: it is now common to fabricate electron device structures with dimensions at 
mesoscopic scale and even at nanoscale. One of the important directions in nanoelectronics 
research is to understand device properties which relate to the existence of small dimensions. 
In this work, we investigate the notion of electrochemical capacitance for conductors in the 
mesoscopic or nanoscale and the nonequilibrium charge distribution at the nonlinear level. 
Using a full first principles quantum theory, a semiclassical theory, as well as a direct numer- 
ical solution, we construct an overall physical picture on the effects of quantum tunneling to 
the nanoscale capacitance. We also investigate the density of state correction to capacitance 
at nonlinear orders of the external bias. For a pair of very small conductors, it has been 
known that the capacitance may be quite different from the usual parallel plate capacitance 
formula given by Co ~ 1/W^, where W is the distance between the two plates. Apart from 
the usual electrostatic fringe effect, there are quantum corrections to the classical formula. 
Quantum corrections come from several sources: a finite density of states of the plates, a 
finite screening length to the electron-electron interactions, and quantum tunneling. 

The quantum correction to classical formula due to density of states (DOS) has been 
theoretically |]I],0] and experimentally investigated in the literature by a number of au- 
thors. For semiconductor hetero junctions they found that DOS contributes a factor to the 
capacitance given by Cdos = e^{dN/dE) where dN/dE is the total density of states of the 
system. Thus the total capacitance C is a result of a series connection of Co and Coos'- 
\jC = 1/Co + 1/Cdos- This formula has been theoretically studied from a dynamic point 
of view and was derived in the low frequency limit of an AC theory [^^ . Significantly, these 
investigations on DOS effects focused on the linear capacitance coefficient C, namely C is 
the linear coefficient of the charge pile-up on a capacitor plate as a function of the external 
bias voltage. Recently some attention [^ has been paid to the nonlinear regime: due to 
the nonlinear bias dependence of local DOS there is also a nonlinear bias dependence of the 



nonequilibriuin charge. The nonhnear capacitance coefficients is one of the topics which will 
be further investigated below. 

Mesoscopic electrochemical capacitance has been found |^ to obey, within a magnetic 
field, weaker Onsager-Casimir symmetry relations. For example it is no longer a symmetric 
function of magnetic field [0 . The asymmetry of electrochemical capacitance has been ob- 
served for a gate close to the edge of a quantum Hall bar |^. The magneto-electrochemical 
capacitance of a three-dimensional quantum dot with three-probes has been studied numer- 
ically in detail in Ref. |^. It is found that at low magnetic fields the magnetocapacitance 
shows a large asymmetry under a magnetic field reversal. At higher fields the capacitance 
is dominated by Aharonov-Bohm type oscillations and the fluctuations of the asymmetry is 
somewhat reduced. For the III-V tunneling heterostructures, the contribution of the density 
of states on the magnetocapacitance is also studied |10[ . The investigation of the frequency 



dependent electrochemical capacitance for a parallel plate capacitor within the nonequilib- 
rium Green's function formalism show interesting oscillatory behavior which is related to 
the retardation effect of the Maxwell equations ||rT| . 

As mentioned above, quantum tunneling effect changes the capacitance value as predicted 
by the classical formula. This effect was recently addressed using numerical analysis of 
atomic junctions ^^. Numerical calculations [l^ of aluminum atomic junctions with tiny 



DOS showed that at small distances W, the electrochemical capacitance C = C{W) actually 
increases with W which is due to tunneling effect. One expects that at larger W when 
tunneling effects is diminished, the capacitance would follow a crossover to the classical 



prediction. However due to the very small DOS of the atomic junction |]T2[, no crossover to 
the classical formula was found in these atomic systems. 

The correction to classical capacitance formula due to a finite screening length was most 
clearly demonstrated from a dynamic point of view on the electrochemical capacitance, 
due to the work of Christen and Biittiker |13| where a conducting quantum point contact 
(QPC) was found to establish a nonequilibrium charge resulting to a finite electrochemical 
capacitance. In particular they have derived a formula for a QPC with a semiclassical 



method [113 



r ~ r ^ p2dNi ^ ^2dN2 ' y'-) 

where R is essentially a reflection probability of the QPC, Co is geometric capacitance, 
dNi/dE and dN2/dE are the total DOS in the regions to the left and to the right of the 
QPC. Qualitatively, the numerical data of the aluminum tunnel junction [Q were consistent 
with Eq. (|l]) in that C is proportional to R. Formula (P is termed "semiclassical" because 
not all the relevant scattering local partial density of states (LPDOS) were included in its 
derivation. The notion of scattering LPDOS was proposed by Biittiker [HI and subsequently 



by Gasparian, Christen and Biittiker |TJ] , and it plays a very important role in low frequency 
AC transport as well as nonlinear DC transport. LPDOS describes the probability of various 
scattering processes [1^. Consider a tunnel barrier as shown in Fig. (1). An example of a 



LPDOS is denoted by da22{Y)/dE which is the contribution of carriers at position r to the 
DOS, and these carriers come from region 2 and ultimately return to region 2. Although 
region 2 is on the right hand side of the tunnel barrier (see Fig. (1)), da22{T^)/dE ^ even 
when r is on the left hand side of the barrier due to tunneling. In deriving |T^ Eq. (|ip for 
a QPC, contributions such as da22{r)/dE with position r on the other side of the QPC, has 
been neglected. 

In this paper, we will further investigate nanoscale capacitors where the two conductors 
have a DC coupling, namely there is a DC "leakage" from one conductor to the other. For the 
linear electrochemical capacitance of a tunnel barrier, we improve formula (|l]) by including 
the tunneling contributions of various LPDOS. This way a full quantum capacitance formula 
is derived and will be compared with (|l]). For a single tunnel barrier there is a quantitative 
difference between these results in the quantum regime, and the difference diminishes as the 
classical limit is approached. The quantum formula and Eq. (0) allow investigations of a 
crossover from tunneling dominated regime to the classical regime, by varying the barrier 
width W. Our derivation as well as the derivation of Eq. (|I]) are within the discrete potential 



model [|l^] that used an approximation where the space is coarse grained into a few regions. 



For the tunnel barrier they are regions to the left of the barrier (denoted by Qi), to the 
right of the barrier (^2), and the barrier region. To confirm that this approximation does 
not affect the predictions qualitatively, we have carried out extensive numerical calculations 
of the LPDOS by directly solving them without the approximation. 

Recently, the theory of non-linear electrochemical capacitance has been formulated using 
the response theory |16|. The electrochemical capacitance of a parallel plate capacitor is a 
nonlinear function of the bias voltage due to the finite DOS near the plates as mentioned 
above. In this work we will study effect of screening on the nonlinear electrochemical ca- 
pacitance for the "leaky" capacitor, which is an important problem not investigated before 
and is relevant for experiments of scanning capacitance microscopy [|1^] applied to nanosys- 



tems. We will derive a general expression of the second and third order nonlinear quantum 



electrochemical capacitance using the discrete potential model |T^. Our analysis naturally 
deduces, in appropriate approximations, a semiclassical expression of the second order non- 
linear electrochemical capacitance for QPC Finally, to compare with results of the discrete 
potential model and semiclassical result, we have directly solved the Poisson equation and 
calculated numerically the linear and the second order nonlinear electrochemical capacitance 
as a function of barrier width of a single tunneling barrier. 

The main results of our investigation are summarized in the following sections. In the 
next section we present our theory of the nonlinear electrochemical capacitance where full 
quantum tunneling effect is taken into account. At the linear order, we compare the quantum 
formula with the semiclassical formula; and using scattering Green's functions we derive 
second and third order nonlinear results. In Sections III we present numerical calculations 
which is compared with the theoretical analysis. Finally the last section summarizes the 
main findings. 



II. THEORY 

In general a two-probe system can be considered as having three regions, a scattering 
region and two electrodes. This is illustrated in Fig. (1) where the scattering region includes 
the scattering potential barrier, and two electrodes are the regions to the left (Qi) and to the 
right (^^2) of the barrier. We are interested in the electrochemical capacitance of this system 
by including the full quantum effects. If we refer regions Qi and ^2 as the two conductors of 
a capacitor, we are dealing with a "leaky" capacitor since the potential barrier provides a DC 
coupling between the conductors. Far away from the the regions, the system is connected 
to contacts which are viewed as large thermodynamic reservoirs, hence in the contacts the 
electron distributions are Fermi-Dirac. When a voltage Vi is applied at contact 1 and V2 at 
contact 2, assume Vi < 0, the electron energy band at contact 1 is changed by dfii = eVi 
and at contact 2 by dfi2 = eV2- The relative electrochemical potential difference is thus 
djji = dni — dn2- due to d^ electrons are injected into the system. The force acting on 
electrons comes from a combination of external and internal fields. In principle, motion of 
electrons in the total field can be solved by Schrodinger equation. In particular we will adopt 
the scattering matrix approach formulated by Landauer |]18[, Imry [jl9[, and Biittiker pO|j2T[] 



to solve the single electron transport problem which gives the necessary LPDOS needed for 
the calculation of electrochemical capacitance. 

Study of electrochemical capacitance is closely related to the calculation of changes of 
the local band ef/(r). It is clear that this local band change near the tunnel barrier is 
different from the shift d^k which occurs at the contacts far away from the barrier. At 
equilibrium conditions the electron energy near the barrier is given by Et = E + d^k — eU{r) 
where E is the electron energy at Fermi level without the applied voltage, dfik denotes 
the electrochemical potential change in reservoir k. Near the barrier electrons accumulate 
for regions where Et > E and deplete for regions where Et < E. It is these accumulated 
charges which we must evaluate. The internal potential build-up ef/(r) can be solved by a 
self-consistent Poisson equation. For simplicity of discussion, in the following we use Ui (r) 
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and U2{r) to denote this potential in regions fli and Q2 respectively. Furthermore, analytical 
derivation of capacitance formula in terms of microscopic quantities is possible if we use a 
space-averaged potential Uk to replace the space dependent potential Uk{r), as was done 



in Ref. [III]. This corresponds to the discrete potential model proposed by Christen and 



Buttiker P|. 



We represent the number of electrons in the region Qk {k = 1,2) incident from contact a 
(a = 1, 2) by aka, which is a function of electron energy Et. Hence cTka = CTkai^E + dua — eUk) ■ 
The number of electrons without external bias (at equilibrium) is thus aka{E), because 
Uk —>■ when rf/i — > 0. By definition, the electrostatic (geometrical) capacitance Co between 
the two regions Qi and ^2 is given by Co = A(5i/(t/i — f/2) (or by Co = AQ2/{U2 — Ui)) where 
AQk {k = 1, 2) is the charge measured from the equilibrium value in region Qk regardless 
where they have come from, i.e. AQk = J2aWka{Et) — aka{E)], where, to avoid confusion 
we use k = I,II to denote the regions from now on. Since there are two electrodes, i.e. 
a = 1,2, AQk thus consists of two parts. For example, in region Qj [i.e. k = I), a. part of 
AQj is due to electrons incident from electrode a = 1 which are scattered back to region 
k = I. We denote this part of AQj by ANi{ni) = an{E + rf/ii - eUi) - (Xn{E). The 
second part of AQj comes from electrons launched at electrode a = 2 but ended up in 
region k = I, this part is expressed by AA^2(^/) = crn^E + d^2 — ef/i) — (Tj2{E). Hence 
AQi = ANiinj) + AN2{nj). 

The above partition of local charge according to where it comes from can be equally 



applied to the scattering local partial density of states [0. Hence, for example, dai2{^i)/dE 
is the LPDOS which is the DOS for an electron incident from electrode 2 passing through 
region Qj and reaching electrode 1. Similarly, da22{^i)/dE is the LPDOS which is the 
DOS for an electron incident from electrode 2 passing through region Qj and eventually 



returning to electrode 2 |2^. Both of these LPDOS describe the tunneling process. This 
latter term is neglected for a semiclassical calculations and is nonzero for a quantum analysis, 
as emphasized in Ref. ^^. They both contribute to the electrochemical capacitance ||2^ 



which is the experimentally measured capacitance defined by 



'^ dfii — dfi2 

The rest of the paper is devoted to calculate C^ including all the quantum effects discussed 
above. 

Based on the above discussions, we can write down the following two equations for 
the classical geometrical capacitance. Using charges of region fij, 

Co X {Ui - U2) = an{E + dfn - Ui) - an{E) 

+ai2{E + dfi2 - Ui) - cTi2{E) . (3) 

Using charges of region Qji, 

Co X (f/2 - f/i) = ajn{E + rf/n - U2) - ajn{E) 

+ain{E + d^2-U2)- cxin{E) . (4) 

Because the same charge defines electrochemical capacitance C^ as given by Eq. (0), we 
have 

Co X {Ui - U2) = C^x {djji - tZ/ia) . (5) 

Finally, it is important to remember that the internal electrostatic potential Uk is a function 
of the electrochemical potential at the reservoirs, 

Ui = UiiHi, H2) , U2 = U2{Hl,H2) ■ (6) 

In above equations we have set electron charge e = 1 so that dfia = Va which is the bias 
voltage at reservoir a. 

Equations (^, |^, ^ are the fundamental equations which we will use to derive quantum 
corrections to Co at the linear and nonlinear orders. Because our theory is gauge invariant, 
without loss of generality we set V\ = V and V2 = throughout the following analysis. 



A. Linear electrochemical capacitance formula 

As discussed above, a semiclassical formula of the linear electrochemical capacitance has 
been derived in Ref. [|13| in the form of Eq. (|l]). In this subsection we derive a full quantum 



formula. 

Taking derivatives of Eqs. (0,^, |5|) with respect to V, we obtain 
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(9) 

where Ei = E + Vi-Ui,E2 = E + V2-Ui,E3 = E + Vi- U2, and E^ = E + V2- U2. In 
deriving the last equation, we have assumed that C^ has no bias voltage dependence |2^ . In 



general the above derivatives should be done at a finite bias voltage V, but experimentally 
one can control this parameter and use very small voltages Q V^ << -E. Hence we will 
evaluate the derivatives at the V ^ limit. In the above equations, the quantity da^a/dEi 
is just the LPDOS in the corresponding regions as discussed above (where we used the 
notation such as daji/dE). 

From Eqs. (|^, ||, H), eliminating dUi/dV and dU2/dV, we obtain 
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(10) 

The electrochemical capacitance C^ calculated from this formula is fully quantum: all the 
tunneling effects are taken into account through the appropriate LPDOS which can be 
evaluated from quantum scattering calculations (see below). 

The general result (p!0| ) can be reduced to the semiclassical form Eq. ([1|) if we apply 
the semiclassical version of the LPDOS. In the semiclassical limit, Ref. |T^ showed that the 
LPDOS is given by 



^ = E^4| + Mi?5,.-|)] (11) 

p 

where T is the transmission coefficient, R is related to the reflection coefficient, Dj = 
daji/dEi + daj2/dE2 and Djj = dam/dEi + da'jj2/dE2 are essentially total local DOS in 
regions Qj and Qu. Substituting (|ll|) into Eq. ([10|) , it is straightforward to prove that 
Eq.(|IOD reduces to the result of Ref. [|13[: 



R _ I 1 1 

where we used notation Cu to denote the /mear electrochemical capacitance C^. If we further 
set R = 1, i.e. no DC coupling is allowed between the two regions, formula (|12|) reduces to 
the familiar electrochemical capacitance of two plates where there is no DC current flowing 
through p. 

In Section III we will provide numerical plots of the full quantum and semiclassical 
formula, and compare them with direct numerical solution of the same problem which does 
not employ the discrete potential model. 

B. Nonlinear electrochemical capacitance formula 

We now derive the second order nonlinear electrochemical capacitance from the funda- 
mental equations (|^), (0), and @. A similar procedure leads to higher order results. To 
obtain nonlinear results we expand Eqs.(^) and (^ order by order in terms of the bias 
voltage Vj3 and internal potential Up. The expansion coefficients are energy derivatives of 
the spectral function cr^a, where the flrst derivative is the linear LPDOS used in the last 
section, while the second derivative is the second order nonlinear LPDOS which is analyzed 
in the Appendix A. Similarly higher order derivatives are the corresponding higher order 
nonlinear LPDOS. To simplify notation in the following we denote Dka = daka/dE and 
Dka = d aka/dE . 

To second order in bias voltage, Eqs.(^ and (|[) become 
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Co{Ui-U2) = Y.Dif5Vp-DiU^ 



/9 

+ Y.\di,{V,-U,Y (13) 



- CoiU, - t/2) = E DiiisVp - DnU2 

+ \T.Diip{Vp-U,f . (14) 

^ /3 

Using Eq.(p!TD and expression (^) of Appendix A, in the semiclassical limit the above two 
equations become 

Co(f/i - U2) = Dj{^ + R)V, + Di^V, - DIU^ 
+ ^RDJ{V^ - U,f + ^TDj{V2 - U^f (15) 

and 

- Co(t/i - f/2) = Djj^Vi + Dn{^ + R)V2 - DjjU2 

+ ^TDii{Vi-U2f + ^RDn{V2-U2f . (16) 

In terms of Cu of Eq. (|T2|), we obtain internal potential f/i and f/2 to first order in 
voltage, 

U, = RV, + ^(Vl + V2) - ^{Vi - V2) (17) 

and 

^2 = i?^2 + ^(^l + ^2) + -^(V^l-V^2) . (18) 

Substituting Eqs. ([l7|) and (|18D into the quadratic terms of Eqs.(|l5]) and ([T6|) , we obtain 



Co(t/i - U2) = Dj{^ + R)V, + Dj'^V2 - DjUi 

+ lTDr{R+'^-^)\V^~V2f (19) 
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and 



Co(f/i - f/2) = Dn^Vi + Dn{^ + R)V2 - DnU2 



1, 



+ -TDij(R+- 

2 ^ 2 D// 



^ ^''?{v,-v,r 



1 



^ C*!! 2/.. .,x2 



+ 2«""( 2 + I^) (''' - ""' 



Combining the above two equations, we finally arrive at 



1 



Co(f/i - f/2) = CniVi - V2) + -Cin(Vl - V2) 



(20) 



(21) 



with the nonlinear capacitance 



Cm — C*ii 
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11 \2 



(22) 



Dn' 2 /},/ 

This result indicates that the second order nonlinear electrochemical capacitance can be 
expressed in terms of microscopic quantities such as the various LPDOS as well as trans- 
mission and reflection coefficients. All of these are calculable and have been studied before. 
Hence this result is very useful in practical predictions of nonlinear capacitance coefficient, 
and it is valid even if there is a DC coupling between the two polarization regions of the 
conductor. 

The general expression (^) is reduced in certain limiting situations. First, for a spatially 
symmetric system where Dj = Djj and Z)/ = Z)//, Eq. (p2D gives Cm = 0. This is expected 
due to symmetry: since Cm is the coefficient of the charge distribution expanded in second 
order of bias voltage {e.g. Eq. (pT])), it must vanish as charge Q -^ —Q for symmetrical 
systems when V — > —V. Second, for a capacitor without DC coupling between the two 
conductors, i.e. for cases T = identically, Eq.(p2D becomes 






(23) 
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which was first derived in a response theory |T6|. Finally, a point worthy some discussion is 



the "resonant transmission point" by setting T = 1 and R = 0. For this case from Eq. (|12|) 
the linear electrochemical capacitance Cn = 0. But from Eq. (521) Cm ^ and is given by 



which is generally nonzero. Apparently we would expect no charge accumulation when T = 1 
hence Cm and all other capacitance coefficients would vanish. However the T = 1 limit 
in the above formula only states the fact that injected charges are going through from one 
capacitor plate to the other at the linear order, and it does not implicate the behavior of the 
charges at nonlinear order where in general T = T{E, U). Thus in setting T{E) = 1 in Eq. 
(p^) is not the true resonant transmission point: at nonlinear order the resonance occurs at 
T{E,U) = 1. 

C. Analysis beyond discrete potential model 

So far we have derived the linear and nonlinear electrochemical capacitance coefficients 
within the discrete potential model, in which the internal potential Uk is parametrized in 
terms of a geometrical capacitance Co- This parametrization is necessary in order to carried 
out analytical derivations, and it is adequate to reveal qualitative features of the physics. 
On the other hand, if one is willing to perform numerical calculations, it is possible to go 
beyond the discrete potential approximation. In this case we can solve the internal potential 
U = U{r) from a self-consistent Poisson equation. In this subsection we derive capacitance 
expressions which are suitable beyond the discrete potential model. 

We start from the charge pile-up written as a three-dimensional spatial integral of the 
charge density [1^ 



Qa= p{,x)d X 

^T.Capyp + lT.CaP,VpV, + ... . (25) 

^ Pi 
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Ref. has shown that charge density p{x) is given in terms of the hnear and nonhnear 
LPDOS, as 






(26) 



2^ dE'' 

To proceed further we must solve the internal Coulomb potential U{x) by the Poisson 
equation 



V^U{x) = Anp{x) 



(27) 



As done previously P,|T6|, for perturbative analysis of the electrochemical capacitance we 
introduce the characteristic potential u{x) 



U{x) =^Ua{x)Va + -J2 UafsVaVp + ... 



(28) 



a/3 



Hence instead of solving U{x) we solve for u{x) order by order. From Eqs.(^) to (p8D, we 
find Poisson like equations for the characteristic potentials [ p6| 



and 



2 da daa 

V M„ + 4:71-— Ua = 4:71—— 

ah ah 



V72 , /I '^^ A ^^"/3 



(29) 



(30) 



where ||T6| 



da. 



a/B 



d'^aa 



al3 



d'^aa d'^ap d^a 



-Up 



-Ur 



(31) 



dE dE^""" dE^-^' dE^'^'' ' rfE^''"''^ 
With the help of Eqs.(|^) and (^), the electrochemical capacitance can be calculated 



from the following expressions, 



Cap = / Qpix)d X 



dan da 



^/3 



Via \ dE dE 



-jr-Up d X 



(32) 



^a/37 = 



iZo- 



{x^d^x 



fdopj 
Qc \ dE 



j^upAdx 



(33) 
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where Qf^^x) and Q/3^(x) are linear and nonlinear nonequilibrium charge distributions. These 
results are useful for numerical calculations where all the quantities on the right hand side 
can be obtained accurately. For instance Eq. (^) has been used in the analysis of atomic 
junctions [|T2|. Eq. (p3| ) is derived for the first time here. 

To end this section we note that in a numerical calculation, the LPDOS daa/dE can be 
calculated using the scattering wavefunction ^7} 



^ - limf (34) 

where v is the velocity of the carrier and ^(x) is the scattering wavefunction for incident 
wave coming from lead a. Eqns. (|2^, |30|) can be numerically solved on a three-dimensional 
grid, for instance a multi-grid technique was employed in Ref. ( [Q) for such a purpose. 

III. NUMERICAL RESULTS 

In this section we present numerical plots for electrochemical capacitance of the tunnel 
barrier in Fig. (1). The numerical curves were obtained along two lines: by plotting the 
analytical expressions (|T0| , |T^ , |2^) which are within the discrete potential model; and by 



direct numerical solution of the self-consistent internal potential U{r) and then applying 
expressions (^) and (|33D . 

To be specific, we choose a numerical calculation box with size x^ — x/j = 12Ai? where Xp 
is Fermi wavelength of the scattering electron. Here xl,r are the positions of left and right 
boundary of the calculation box. We fix the tunnel barrier of width W symmetrically in the 
center of the calculation box. This way the regions Qj and Qu discussed above are given 
by the space between the calculation box and the barrier walls. The quantum scattering 
problem by this single barrier is easily solved, from which we obtain various LPDOS using 



the scattering wavefunction according to Eq. p^) . To apply expressions (|iy, |T2|, |22|) which 
are appropriate for the discrete potential model, we spatially average these LPDOS in the 
respective regions which gives us the corresponding quantities in these expressions. On the 
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other hand, in applying expressions ( p2| , |33[) which uses the full spatial dependent internal 
potential, the spatial integration range should be the Thomas-Fermi screening length [Q as 
discussed in Appendix B. The screening length is determined |^ by solving the Poisson-like 
equations ( p9D (and (|30| ) in the nonlinear case). From now on we will use atomic units such 
that h = 2m = e = 1. In typical nanoscale systems P3 with charge density 10^^, Fermi 



wavelength A^ ~ 47nm. In the following we use A^ to set the unit for length and choose 
Fermi energy Ep = 0.31. 

Fig. (2) plots the linear electrochemical capacitance obtained from different approaches 
as a function of the barrier width W for the fixed barrier height Hq = 0.8: (a), the electro- 
chemical capacitance C calculated numerically from Eq.(|32|) (solid line); (b). the analytic 
quantum electrochemical capacitance formula in the discrete potential approximation Cg 
from Eq.([TO|) (dotted line); (c). the semiclassical electrochemical capacitance [jl3[ Cg of Eq. 



(|1^) (dashed line); (d). the effective classical geometric capacitance Co ~ 1/W (dot-dashed 
line). For very large the barrier width, it is clearly shown that all curves approach to the 
classical behavior ~ 1/W since quantum tunneling effect is negligible. For thin barriers 
where tunneling effect is significant, the behavior of electrochemical capacitances C, Cg, 
and Cg are completely different from the classical regime. In this quantum regime as one 
increases the barrier width, the electrochemical capacitance increases rather than decreases. 
This increasing behavior at very small W is expected since tunneling tends to diminish 
charge polarization, thus C ~ when VT ~ 0. Hence C{W) should indeed start from small 
values and increase a bit before it goes down when W is large enough. 

To examine the DOS correction which is another quantum effect, we note that one can 
only separate out the geometrical effect from the DOS effect in the semiclassical limit (as in 
Eq. (p!2D), and in general these effects are mixed. Furthermore, in a discrete potential model 
all the quantities (both in quantum and semiclassical calculations) are spatially averaged, 
hence capacitances are under-estimated. This is why both Cg and Cs curves are consistently 
below the full numerical solution C. Fig. (2) shows some difference between the quantum re- 
sult Cg and semiclassical result Cg- To understand this difference we have plotted the partial 
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DOS dnii{Qjj)/dE (solid line) and dni2i^ii)/dE (dotted line) in the inset of Fig. (2). As 
expected, dnii{Qjj)/dE goes to zero for large barrier widthes where the semiclassical theory 
is a good approximation. It is nonzero in the quantum tunneling regime for small barrier 
width. dnii{Qjj)/dE is also numerically much less than dni2{^n)/dE. Hence neglecting 
dnii{Qji)/dE in the semiclassical analysis gives a small difference between Cs and Cg in the 
tunneling regime (see Fig. (2)). To further compare with the semiclassical result of QPC of 



Ref. [|13[, we have also examined the behavior of capacitance by varying the barrier height 
Ho for a fixed barrier width W: the results using Eqs. (^,|l^,|l^) are, again, similar in the 
quantum regime and the same in the classical regime. When the barrier height Hq is rela- 
tively small, the appearance of quantum mechanism leads to a correction for semiclassical 
electrochemical capacitance. 

The physical behavior of second order nonlinear electrochemical capacitance coefficient 
Cm can be studied for an asymmetric barrier: as discussed above Cm = for symmetric 
systems (see Eq.(^2])). For the asymmetric barrier where the barrier heights are Hi = 0.2 
and H2 = 1.0, shown in the inset of Fig. (3a), Fig. (3a) plots Cu versus W and Fig. (3b) plots 
Cm. The linear coefficient Cu is very similar to that of Fig. (2) of a symmetric barrier. For 
the full quantum numerical result (solid line). Cm starts from zero and becomes negative 
for very thin barrier, reaches minimum at around W = 1.0, and then oscillates around zero. 
The oscillatory behavior can be traced to oscillations in second order DOS D = d'^N/dE'^ 
of Eq. (p2D . In the inset of Fig. (3b), we plot PDOS Dj and Dj. Correlating the behavior of 
Cm and PDOS, it is clear that the negative values of Cm is due to D. 

In Fig. (4), the linear and nonlinear nonequilibrium charge distribution for this asymmet- 
ric barrier, Qi{x) and Qn{x), are plotted. These quantities, especially Qii{x), have not been 
studied carefully before. It is thus interesting to offer several observations, (a) The linear 



charge distribution Qi{x) is in the form of a resistance dipole |^, whereas the nonlinear 
charge Qu{x) is more like a quadrupole. (b). The linear charge distribution is numerically 
much larger than the nonlinear charge distribution. The total charges are conserved, i.e., 
J Qi{x)dx = J Qii{x)dx = 0. (c). In the discrete potential model, the average nonlinear 
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charge Qu is numerically even smaller. Because of this spatial average, the nonlinear charge 
distribution becomes a dipole in the discrete potential model. This is responsible for the 
difference between full quantum calculation and that of the discrete potential model. 

IV. DISCUSSION AND SUMMARY 

In this work we have investigated the quantum version of a "leaky capacitor" in the 
coherent nanoscale regime in both linear and nonlinear order in terms of the external bias 
voltage. We have derived an analytic formula of electrochemical capacitance where the two 
plates have a DC couphng, and tunneling effects between the two plates are fully included 
by explicitly using partial local density of states duaai^k) / dE . Within the discrete potential 
model where all quantities are averaged over the polarization regions, analytic expressions 
for the linear and second order nonlinear electrochemical capacitance have been derived. In 
addition, linear and nonlinear capacitance coefficients formula are derived in terms of the 
self-consistent potential, and these formula are suitable for numerical calculations in the 
full quantum regime. Our calculation showed that the analytic results using the discrete 
potential model are consistent with the full numerical solution, for the single tunnel barrier 
structure. The electrochemical capacitance formula derived in this paper are suitable for 
analyzing AC transport at relatively low frequency. At very high frequency, one has to 
consider radiation effect thus the quantum equation must be solved self-consistently with 
the full Maxwell equation instead of the Poisson equation used here. 

Quantum behavior of the electrochemical capacitance is found in the tunneling regime 
that the capacitance increases with the barrier width W. This is in clear contrast to the 
classical behavior of 1/W. What is the condition that this non-classical phenomenon be 
observed ? Let's consider this question using the semi classical formula [|I^ Eq. (|1^) which 
can be rewritten as 

Co ~^ D 



For tunneling, R ^ [1 — exp{—W/l)] where / is a characteristic length depending on system 
details such as the barrier heights. When C = C{W) increases with W, we have dC/dW > 
which gives to a condition for the range of W that allows the nonclassical behavior. Using 
Eq. (p5D, for tiny capacitor plates D « Co-, one can have a reasonable and experimentally 
accessible range of order I. On other hand for large plates D » Co, one can not observe 
the nonclassical effect unless W is several orders smaller than / which is not experimentally 
accessible. Hence the condition to observe the non-classical behavior is tunneling and also 
small DOS. Systems which satisfy these conditions are nanoscale capacitors, whereas capac- 
itors with large plates such as Josephson junctions (macroscopic) do not satisfy the DOS 
condition. 

The nonlinear theory developed here can be pushed to higher order analytically within 
the discrete potential model. At linear order the full quantum formula ([lOD and the semi- 
classical formula (|T^) give certain numerical difference in the quantum regime. There is a 
more substantial numerical difference between the discrete potential model and the full self- 
consistent numerical calculation using expressions ( |5^ and (P5|), although all these results 
are qualitatively consistent. At second nonlinear order, the nonequilibrium charge distribu- 
tion behaves as a quadrupole, this is compared to the resistance dipole when linear order 
charge is considered. It is interesting to note that the nonlinear charge can be nonzero when 
linear charge is zero: this happens at the linear resonance point. The nonlinear capacitance 
coefficient also behaves quite differently from the its linear counterpart, as shown by its 
oscillatory behavior linked to the nonlinear LPDOS. 
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V. APPENDIX A. 

The nonlinear electrochemical capacitance depends on the nonlinear LPDOS, as shown 
in Section II. In this Appendix, we derive the nonlinear (2nd order) LPDOS using Green's 
functions. In particular we have to relate the second order LPDOS dPaap/dE'^ to the total 
LDOS (Pa/dE"^, where indices a, j3 label the leads. 



We start from the definition of LPDOS expressed in terms of the Green's function |30 

^^^ = Re[2miG'r2G''T,G%.] (36) 

where G^ is the retarded Green's function, Pq, is the linewidth function, and i?e[...] denotes 
the real part of [...]. Using relation |M 



we have 



{G''MG'')xx — 2^ Glxi^xiXiGl^x 

X\X2 

X\X2 

= Gl,Tr[G^M] (38) 



where M is a matrix. Eq. (|36[) becomes, 

^ = -2nIm[G:,Tr{r,G'^r^Gn] 

= ^(GL-GL)T = -— (39) 

where T = Tr(P2G'"riG"')/(47r^) is the transmission coefficient |2^ which is a real quantity; 
da/dE is the total local DOS. 

Taking energy derivative of Eq. (p6D , we obtain 

^ = 2TTlm[iG'-G'T2G''T,G%x 
+ (GTsG^G^PiGOx. 
+ (GTsG^PiG^'G^....] . (40) 
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The first term of (^OD can be simplified using Eq.(p8D as follows, 

= GlJr{GnTr{V^G-V,G'-) = l^^G^G^, . (41) 

The other two terms of Eq. ( |1(]| ) can be simplified in a similar fashion. We thus have 

^ = I-Im[2{G^G'U + (G^G")..] 

= -^[{G'G^U - [G'^G'^U] = T^ . (42) 



In deriving the last equation we used the fact that G^G" is a real quantity. Using (|4^) we 
find 

(i^cTii d'^(J22 (Per d'^CTu d?<y2i 
dE^ ^ dE^ dE^~ dE^ ^ dE^ 

-i^-^T)^- (43) 

Now we consider a system with a DC current passing through. Due to polarization 
we again divide the system into two regions fi/ and Vtn. In the semi classical treatment 
where the tunneling is neglected, the partial DOS d'^a22/dE'^ = in region I, and similarly 
d'^aii/dE'^ = in region II. These relations and Eqs.(^^ and (^31) lead to 



%^ = ^AT + ^./^[(l - 2T)S^, - T]} (44) 

where k labels the polarization region Qk- For two probe system, it gives 

^^§|r+(i-2r)A„j, (45) 

This expression is the semiclassical second order nonlinear LPDOS, which is in contrast to 



the semiclassical linear LPDOS Eq.(|TT]) derived in Ref. |T3[. The nonlinear LPDOS plays a 
crucial role in determining the nonlinear electrochemical capacitance, as given in Section II. 
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VI. APPENDIX B 

In this appendix, we give an example of calculating the second order nonlinear capaci- 
tance Cm by directly solving Poisson equation. This can be done analytically only for very 
simple systems. 

Consider a system which consists of three regions: two leads (region I and III) and 
an infinite potential barrier (region II). Since the calculation is perturbative, we have to 
calculate the linear characteristic potential by solving Poisson equation Eg . (|29|) . We assume 
that the partial local DOS dai/dE and (Pai/dE"^ are constant in region I and zero otherwise 
^. Similarly da2/dE and d'^a2/dE^ are constant in region III and zero otherwise. The 
solution of the Poisson equation Eq. (^) is 

X 

region I : Ui = 1 — Aiexp{—) 
region 11 : Ui = aix + hi 

X 

region III : ui = y42exp(— — ) (46) 

A2 



where Ai, A2, ai, and 61 are constants to be determined. In Eq.(^6D, we have defined the 
screening length A~^ = Andaa/dE and the boundary conditions that -ui ^ 1 as x -^ — cxo 
and Ml — i> as X -^ 00 have been used. Using the boundary condition that ui and dui/dx 
be continuous at x = a/2 and —a/2, it is straightforward to find 

1 ai 

Qi = I X I X > 61 = — a + 2A2 
a + Ai + A2 2 

Ai = aiAiexp(— — ) , A2 = aiA2exp(— — ) (47) 

ZAi ZA2 

The linear electrochemical capacitance can be obtained immediately, 

dp{x) 



_ f Op{x, 



■dx 



— / V'^Uidx 

A 1 

47r a + Ai + A2 
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(48) 



where A is the cross-section area of the metalhc wire. Using the global DOS dNa/dE = 
Volume daa/dE = XaAda^/dE = A/inXa, we arrive at the result first obtained by Biittiker 

1 Ana 11 

Chi^~^ dNi/dE ^ dN^/dE ^ ^^ 

With the solution of mi, the Eq.(^) becomes 



region I : 


-^1 


1 .2 ,2x 
-2Aexp(^ 


region II : 


V^un = 




xgion III : 


^2 1 
V Mil + ,2^11 = 
^2 


1 .2 /2x 



(50) 

where we have introduced another screening length A~^ = And'^aa/dE'^ corresponding to 
LPDOS d^aa/dE'^ and Ai and A2 are known from the calculation of ui. The solution of 
Eq.(PD]) is 



region I : Uu = -^^? exp(— ) + Bi exp(— ) 



^A?exp(— ) + 5iexp(- 

region 11 : un = a2X + 62 

A 2x X 

region III : Uu = -^A2exp(— ) + 52exp(— ) (51) 

0A2 A2 A2 



After matching boundary conditions at x = a/2, —a/2, we obtain 



The second order nonlinear electrochemical capacitance Cm is 

d'^p{x) 

h 
-1, 






_ 1 /■ 9^p(: 
^''' = 2k^Vi 

-VMll|_a/2 ■ A 



A-n 

Ai 2 / '^ \ -^1 / 
rAiexp(-— ) - — --exp( 



67rAf ' "' \i AtxXi "' 2Ai- 
A 1 .Af A^ 

47r(a + Ai + A2)3^3A? 3Ai^ 

From the definition of the screening length, we have 
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(53) 



A. 



1 



(54) 



A AnAXadaa/dE dNa/dE 

where we have used the fact that there is charge polarization only in the region AX^- Simi 
larly, we obtain 



A. 



A d^No./dE'^ 
With the help of Eqs.([5^), (|^), and (|49|), we finally have 

d'^Ni/dE^ d'^N2/dE^ 



Cm - — 



{dNi/dEf {dN2/dEf 



which agrees with Eq.(p3|). 



(55) 



(56) 
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FIGURE CAPTIONS 

Fig. (1) The energy band configuration for a model barrier. Near the barrier the band bottom 
is different from that away the barrier. The inset is an ideal experimental setup of 
which the energy band is schematically shown in the figure. 

Fig. (2) The linear electrochemical capacitance as a function of barrier width W for a square 
barrier with the barrier height Hq = 0.8. The solid line is the full quantum numerical 
calculation C, the dashed line and dotted line are from the quantum result Cg and 
the semi-classical result Cg in the discrete potential approximation, respectively. The 
dot-dashed line is the classical result C ~ 1/W. The inset: the corresponding partial 
density of states versus the barrier width W. The solid line is dNii{Qjj)/dE and the 
dotted line is dNi2{VLu)/dE. 

Fig. (3) (a). The linear electrochemical capacitance as a function of barrier width W for the 
asymmetric barrier (see inset). The system parameters are Wi = W2, Hi = 0.2, 
H2 = 1.0. (b). The second order nonlinear electrochemical capacitance versus W. In 
Fig. (3a) and (3b), the solid, dashed, and dotted lines correspond to C, Cg, and Cg, 
respectively. In Fig. (3b), we have multiplied Cq and Cs by a factor of 5. The inset 
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of Fig. (3b): the corresponding partial DOS Di (solid line) and Di (dotted line) as a 
function of W . For illustrating purpose, we have divided Dj by a factor of 30. 

Fig. (4) The linear (solid line) and nonlinear (dashed line) charge distribution for the asym- 
metric barrier, where W = 20, Hi = 0.2, and H2 = 1.0. The dotted line shows the 
shape of asymmetric barrier. 
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